Broken Symmetry in Density-Functional Theory: Analysis and Cure 
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We present a detailed analysis of the broken-symmetry mean-field solutions using a four-electron 
rectangular quantum dot as a model system. Comparisons of the density-functional theory predic- 
tions with the exact ones show that the symmetry breaking results from the single-configuration 
wave function used in the mean-field approach. As a general cure we present a scheme that sys- 
tematically incorporates several configurations into the density-functional theory and restores the 
symmetry. This cure is easily applicable to any density-functional approach. 
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The nanoscale semiconductor systems are technically 
very promising for future components of microelectronic 
devices. From the theoretical point of view, quantum 
dot (QD) systems are valuable source of novel quan- 
tum effects. Many of these result from the fact that 
the electron-electron interaction and external magnetic 
field have greatly enhanced effects compared to atoms 
and molecules. This raises new challenges for the theo- 
retical methods, and the validity of approximations in, 
e.g., mean- field approaches can be questioned. For this 
reason, QD systems serve as perfect test cases to develop 
the electronic structure methods, with the results still 
applicable to great variety of physical problems where 
mean-field approaches have been used. 

In earlier studies, Hartree-Fock (HF) and especially 
density-functional theory (DFT) methods have shown to 
produce accurate results for various QD systems, even 
with small N. However, in the context of solutions with 
a broken spin symmetry, the validity of the mean-field ap- 
proaches has been actively discussed in the literature [2j . 
The spin-density wave (SDW) formation in QDs has been 
compared to similar phenomena found in isotropic met- 
als 0, organic linear-chain compounds atomic nu- 
clei 6], and small fermion clusters 0. According to the 
Jahn- Teller theorem, any non-linear molecular system in 
a degenerate electronic state becomes more stable by re- 
moving the degeneracy and thus lowering the symmetry 
and the total energy. A crucial difference between molec- 
ular and QD systems is, however, that as the nuclei in 
molecules are free to move and relax, the QD potential is 
external and fixed as it results from, e.g., metallic gates. 
Thus to lower the symmetry in QD, the spin densities 
must "relax" in an anti-ferromagnetic fashion to a SDW 
solution. This is claimed to reveal the electron correla- 
tions inherent in the true ground state 0). 

In this Letter, we analyze symmetry breaking in a two- 
dimensional rectangular QD 1] using both DFT and ex- 
act diagonalization (ED). We concentrate on the four- 
electron case, as it is the first particle number showing 



the general features of electronic structure seen also for 
larger particle numbers, such as the transitions between 
the two spin states S — and I and the SDW solution 
predicted by DFT. We find that SDW clearly reflects the 
limitations of basic DFT to describe systems that have 
more than one major configuration in the ground-state 
wave function. There is a continuous interest for devel- 
oping DFT methods for this kind of systems. The main 
difficulty for DFT is the fact that these systems have 
ensemble- v-representable densities (E-VR) in contrast to 
the more common pure-state w-representable densities 
(P-VR) 3J. As an interesting feature we see a contin- 
uous transition from an E-VR to a P-VR density as we 
deform our QD. Finally, we present a simple modification 
of DFT that is able to describe the multi-configurational 
nature of the ground states. 

The generally used model Hamiltonian of an N- 
electron QD system can be written as 



H 



N 

E 



2m 



■Vf + WrO 



N 

E 



er, 



(1) 



where we have used the effective-mass approximation to 
describe electrons moving in the xy plane, surrounded by 
background material of GaAs with the effective electron 
mass to* = 0.067m e and dielectric constant e = 12.4. 
We use scaled atomic units, and energies are thus given 
in Ha* « 11.86 meV and lengths in a* B « 9.79 nm. The 
external confinement in the xy plane is described by an 
infinite hard-wall potential, 
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0, < x < 0L, < y < L . , 
oo, elsewhere. 



The deformation parameter (3 defines the ratio be- 
tween the side lengths of the rectangle. The area of 
the dot is fixed to be it 2 . The single-particle eigen- 
states are sine- functions in both directions, labeled with 
two quantum numbers (n x ,n y ), and energies E n:c ^ ny = 
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{n%./ (3 + P n y) /2. Fig. shows the three lowest eigen- 
states and the most important S z = configurations of 
the four-electron QD. 




FIG. 1: Left panel: The three lowest single-particle states and 
their quantum numbers (n x ,n y ). Right: Electron occupations 
for the four important S z = configurations C;. 

We solve the electronic structure of QD using ED and 
DFT. In ED, the many-particle wave function is con- 
structed as an expansion of the non-interacting eigen- 
states. The results approach the exact ones as more 
terms are added to the expansion. We use a basis of 
up to around fifteen thousand configurations. The inter- 
action matrix elements are calculated numerically using 
Gaussian integration. In ED, all many-body quantum 
effects are taken into account in an exact fashion. In 
DFT, these are incorporated in a mean-field fashion as 
an effective potential. In the DFT method used, we allow 
different spin densities for up and down electrons. This 
is necessary for 5^0, and needed also for 5 = in order 
to find broken-symmetry solutions. More details of the 
DFT method and the numerical implementation can be 
found from Ref. Q] and references therein. 

In Fig. |21 we present the DFT and ED energies of the 
rectangular quantum dot as a function of the deforma- 
tion parameter /3. For (3 close to unity, the 5 = 1 state is 
lower in energy than the 5 = state, in accordance with 
Hund's rule. In the case of the 5 = 1 state, the DFT 
energies compare quite well with those obtained by ED: 
the deviation between the two remains nearly constant 
for all values of (3. Such a behavior is not seen in the 
5 = results, for which we show two DFT energies: one 
with retained symmetry and another with a broken one. 
The broken-symmetry solution has a non-zero total spin 
density, corresponding to a SDW solution, see Fig|3] The 
DFT calculation with the spin symmetry does not con- 
verge for smaller (3 than those shown. This is due to the 
degeneracy in the system. Convergence can be achieved 
by use of fractional occupations. Comparing the DFT 
5 = energies to the exact ones, one can see that, unlike 
for the 5=1 case, the error in DFT is not constant. 
The energy of the symmetry-restricted state grows lin- 
early towards /? — ► 1 where the ED energy saturates. 
On the other hand, the SDW state has an energy that 
overcompensates the error in the symmetry-restricted en- 
ergy. The energy of the SDW state is closer to the exact 
value than the energy of the proper-symmetry state. One 
should note that the errors in DFT energies nearly can- 
cel at the ground-state transition point, and the DFT 
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FIG. 2: Energy of the four-electron dot as a function of the 
axis ratio j3. The solid lines present ED energies, we use 
crosses for 5=1 and diamonds for 5 = 0, correspondingly. 
The dashed lines are DFT energies, pluses for 5=1, boxes 
for the S = SDW solution, and circles for the symmetric 
S — energy. 



prediction for it is very accurate. 

It is claimed that the SDW spin densities reflect the in- 
ternal structure of the system |2| . To analyze this claim, 
we have plotted the SDW spin density of the DFT and 
ED conditional densities in Fig. [3J The conditional den- 
sity is defined to be the electron density of the remaining 
three electrons as the coordinates of one of the electrons 
are fixed. In addition, we plot the ED spin density for the 
sum of the 5 = and 5 = 1 states. One can see from the 
densities that there is a clear anti-ferromagnetic order in 
the system. Densities for parallel spins are localized in 
the opposite corners. Apart from this fact, the similarity 
of the conditional densities to the SDW density of Fig. El 
is marginal. However, the similarity of the SDW den- 
sity to the unphysical mixture of the two spin states is 
very clear. The only difference is that the DFT density 
is slightly more localized. One should note that this sim- 
ilarity of the SDW solution to a mixture of two different 
spin states is pointed out by Hirose and Wingreen using 
ED in restricted basis Q . 

To understand the electronic structure of the system 
and to analyze the problem associated with the SDW so- 
lution, it is enough to consider only the most important 
configurations in the ED solution, presented in Fig ^ 
The 5=1 state, the ground state for small /3, has three 
different S z states which are degenerate in energy. The 
one with S z = consists of configurations C3 and C4 with 
equal weights for all values of (3. The 5 = state is the 
ground state for large j3, and it consists of the configu- 
rations C\ and C2. For (3=1, these have equal weights, 
but for larger (3, C2 moves higher in energy and has a 
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smaller weight in the exact wave function. For (3 « 1.2, 
Ci is clearly the dominating configuration. One should 
note that at this value of (3, proper symmetry is restored 
in the DFT solution. The most natural reason for the 
occurrence of the SDW solution is that the basic DFT 
is unable to take into account more than one important 
configuration for the construction of the Kohn-Sham or- 
bitals and the resulting densities. In terms of the configu- 
rations, the DFT spin densities at (3 = 1.2 correspond to 
C\. For smaller j3, however, the SDW spin densities can 
only be obtained by a linear combination of all four con- 
figurations. For (3=1, this linear combination is equal to 
the unphysical mixture of two different spin states used 
for Fig. |2| (b) above. 

It is possible to analyze the broken-symmetry solution 
more generally by considering a mean-field-type single- 
configuration wave function for two up and two down- 
spin electrons, occupying the orbitals tpa and sm{9 a )ij}\ + 
cos(0 cr )^2, where 6 a contains the variational freedom for 
a spin type a. Expanding this wave function results in 
four configurations similar to {Ci}f =1 above. Assuming 
a further similarity to the QD case for j3 = 1, one can 
write a Hamiltonian matrix of the four configurations as 




FIG. 3: (a): DFT spin density for the SDW solution. The 
density for the other spin type can be found through rotation 
by 90 degrees, (b): ED spin density for the artificial sum of 
S = and 5 = 1 states, (c): Conditional density from ED 
for the same spin type as the electron fixed at '+'. (d): ED 
conditional density for opposite spins. The number of contour 
lines (drawn at uniform spacing) is fixed to 10 in each figure 
to ease comparisons. The SDW density is more similar to the 
unphysical ED density of (b) than the conditional densities. 
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where the configurations couple via the off-diagonal ma- 
trix element S (taken to be real) . The four exact energies 
are Eq ± S and E\ ± S. One can set without loss of gen- 
erality Eq = and E\ = 1. The single-configuration 
energies have an interesting dependence on <5, shown in 
Fig. We present the energy as a function of the two 
variational angles 9 for cases 6 = 0.2 and 0.8. For small 
S, the second orbital for the minimum-energy solution is 
tpi for one spin type and ip2 for the other. For the case 
5 = 0.8, the minima are found with orbitals ipi + ip2 and 
ipi — V'2 ■ The resulting total wave function of this case 
can easily be found to be a sum of the two exact wave 
functions (with energies —S and 1 — S), and the energy 
of the mean-field state is equal to the average of the two 
exact energies. Furthermore, if one assumes that ip\ has 
a node on the x axis and ip2 on the y one, one can find 
densities similar to the SDW solution above. Now for 
the QD, the value of 5 is close to 0.8, and one can un- 
derstand the occurrence of the SDW solution from this 
more general argument. 
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FIG. 4: Single-configuration energy as a function of the two 
angles in the wave function. The left panel corresponds to 
S = 0.2, and the right one S = 0.8. Black areas are the 
lowest in energy. The proper symmetry of the wave function 
is found on the dashed diagonal line. The broken-symmetry 
energy minima of S = 0.2 correspond to a single configuration, 
and for 5 = 0.8 to SDW solutions. 

The SDW solutions of Ref. for larger particle num- 
bers can equally well be understood based on the four- 
electron case and the general argument presented above. 
In the cases where two Kohn-Sham orbitals are degen- 
erate, we have a S = 1 ground state. When the aspect 
ratio (3 is changed, the energies split and one always finds 
a broken-symmetry SDW solution. The similarity to the 
four-electron case follows from the fact that in all these 
cases, there are two spatial orbitals of both spin type oc- 
cupied by two electrons. The spin density in the SDW 
structure can then be directly found from the two nearly 
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degenerate states as (ipi +1P2) 2 — (j^i — V^) 2 - For example, 
the densities in Fig. 6 of Ref. [l| are accurately repro- 
duced by this formula using for the degenerate states the 
non-interacting ones with quantum numbers (1,3) and 
(3, 2) for the left panel, or (3, 2) and (4, 1) for the right 
panel. 

Based on the results presented above, it is clear that 
standard DFT is not able to describe accurately E-VR 
systems. The method of Ullrich and Kohn Q is one pos- 
sible solution, but this method might have an underlying 
problem. Namely, even in the case of an open shell and 
degeneracy, there are systems that still are P-VR, simply 
because the configurations do not necessarily mix even if 
they are degenerate. One such example is the parabolic 
QD, where the angular momentum is a good quantum 
number and single-particle states can be chosen in such 
a fashion that only one major configuration is found. It 
is not straightforward to see how the method of Ref. Q 
assort the open-shell cases that are E-VR from those that 
are still P-VR. 

As a possible solution we propose a scheme where first 
a standard DFT calculation is performed for the system 
(without symmetry breaking but with fractional occupa- 
tions for the degenerate levels) . For our case with (3=1, 
the occupations of the two highest orbitals ipi and ^2 are 
1/2 and the DFT energy is 13.26. One can construct two 
S = S z = DFT configurations that have density equal 
to the DFT one by defining new orbitals <p± = iji-y ± iip2- 
Now the configurations involve the core DFT orbital ipo, 
and either <f> + or 0_ . The occupied orbitals in config- 
urations are the same for both spin electrons, similarly 
to d and C 2 in Fig. □ The coupling of these two DFT 
configurations can be approximated by 

5= [ <P*+(r-x)r+(^) — Mri)Mr 2 ) dr 1( lr 2 . (4) 

J 7*12 

Now the DFT energy gives the diagonal Hamiltonian 
matrix elements and 5 the off-diagonal ones. The two- 
ensemble DFT energy can be found by diagonalizing the 
Hamiltonian matrix. For our example, the absolute value 
of delta is found to be rj 0.14. Thus the mixing of the 
two configurations lowers the DFT energy to 13.12. This 
value is consistent with the ED one, if one takes into ac- 
count the difference in DFT and ED energies for the non- 
degenerate cases. One should note that for a parabolic 
QD, the absolute value of 5 is zero (resulting from the 
rotational symmetry), meaning that the configurations 
(with different angular momentum) do not mix. This 
also shows that our scheme correctly predicts the sys- 
tem to be P-VR. A similar behavior can be found for the 
multiplets in open-shell atoms. 

The underlying idea of the scheme presented above is 
that DFT is able to efficiently describe correlation effects 
in a certain subspace of the full Hilbert space. This sub- 
space is related to a one DFT configuration. In P-VR 



cases this is sufficient for the accurate description of the 
system, but for an E-VR case, there are two or more sub- 
spaces relevant for the ground state, and DFT is unable 
to couple these. This coupling can be introduced, and 
one natural way is via S of Eq. Q above. 

The generalization of the scheme for cases without an 
exact degeneracy of the DFT orbitals is straightforward. 
In addition, the approximation made for S can be directly 
used for cases with larger particle numbers, too. This is 
because the states that are occupied in both configura- 
tions do not appear in the formula for 5. We believe 
that the presented approach shows to be useful for many 
applications of DFT, especially for molecules, where the 
calculations of chemical reactions have observed similar 
problems of basic DFT [°J. More details and compar- 
isons with other ensemble DFT approaches are left for 
forthcoming studies. 

Concluding, have shown that the use of a single- 
configuration wave function in a mean-field theory can 
lead to an unphysical solution with a broken symmetry. 
In our case of a four-electron rectangular QD, the energy 
of the SDW solution is reasonable, but the spin densities 
have only a minor similarity with the exact total or con- 
ditional ones. We also present an analysis with a more 
general Hamiltonian matrix and we feel that our findings 
are relevant for a great variety of systems studied by the 
mean-field approaches, DFT in particular. As a cure, we 
propose a scheme for incorporating systematically several 
configurations into a mean-field approach. The method 
presented avoids the necessity of symmetry breaking, and 
has a built-in criteria to determine if several configura- 
tions are actually needed or not. 
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